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We present new limits on Ve(V e ) — > v T (p T ) and v e (pe) — ♦ Vs oscillations by searching for v e disap- 
pearance in the high-energy wide-band CCFR neutrino beam. Sensitivity to v T appearance comes 
from t decay modes in which a large fraction of the energy deposited is electromagnetic. The beam 
is composed primarily of v tL {y IJ \ but this analysis uses the 2.3% Ve(V e ) component of the beam. 
Electron neutrino energies range from 30 to 600 GeV and flight lengths vary from 0.9 km to 1.4 km. 
This limit improves the sensitivity of existing limits for v £ — > v T at high Am 2 and obtains a lowest 
90% confidence upper limit in sin 2 2a of 9.9 x 10 -2 at Am 2 ~ 125 eV 2 . 

PACS numbers: 14.60.Pq, 13.15.+g 



Neutrino oscillations occur if neutrinos have non-zero 
mass and mixing. If recent evidence Jl],^| for neutrino 
oscillations is confirmed, it will radically alter our under- 
standing of both particle physics and cosmology. Neu- 
trino oscillations may also explain the observed deficit 
of neutrinos from the sun. In the two-generation mixing 
formalism, the oscillation probability is given by 



sin 2a sin 



1.27Am 2 L 



(1) 



where Am 2 is the mass squared difference of the mass 
eigenstates in eV 2 , a is the mixing angle, E v is the in- 
coming neutrino energy in GeV, and L is the distance 
between the point of creation and detection in km. 

While the high Am 2 regions of parameter space for 
v„ — > v P and v,, — > v T oscillations have been excluded 



to very low mixing angles (2 



10 3 for v u 



and 



5 x 10 3 for —* v T ) the high Am 2 



v T parameter 



space is much less constrained as a result of the difficulty 
in producing high energy v e beams. The CCFR sample of 
over 20,000 charged-current u e interactions comprises the 
largest sample of high energy v e 's to date. Previous high- 
statistics limits were obtained from reactor experiments 
with much lower beam energies. Accelerator limits 
were obtained by BEBC j§ and Fermilab E531 § which 
searched for v T appearance in emulsion. 

We previously reported a limit on — > v e oscillations 
by searching for v e appearance in the v e N charged- 
current data sample. In this report we present new limits 
using the same data sample on v e ~ > v r and v e —* v s 



oscillations. Both limits use a v e disappearance test. The 
v e — > v T limit is also sensitive to v r appearance through 
t decay modes in which a large fraction of the energy 
deposited is electromagnetic. 

The CCFR detector [|g consists of an 18 m long, 
690 ton target calorimeter with a mean density of 
4.2 g/cm 3 , followed by an iron toroidal spectrome- 
ter. The target consists of 168 steel plates, each 
3m x 3m x 5.15cm, instrumented with liquid scintillator 
counters placed every two steel plates and drift chambers 
spaced every four plates. The separation between scin- 
tillation counters corresponds to 6 radiation lengths, and 
the ratio of electromagnetic to hadronic response of the 
calorimeter is 1.05. The toroid spectrometer is not di- 
rectly used in this analysis which is based on the shower 
profiles in the target-calorimeter. 

The Fermilab Tevatron Quadrupole Triplet neutrino 
beam is created by decays of pions and kaons produced 
when 800 GeV protons hit a production target 1.4 km 
upstream of the neutrino detector. The resulting neu- 



trino energy spectra for Vu 



and v P are shown in 



Figure |l|. The 2.3% v e component of the beam used in 

this analysis is produced mainly from K ± — > 7r°e ± Vp 
occuring in the 0.5 km decay region just downstream of 
the production target. The v r content of the beam is less 
than 10- 5 . 

Neutrino interactions observed in the detector can be 
divided into three classes depending on the type of in- 
coming neutrino and interaction: 



1. v M N 



X (y^ charged-current (CC) events). 
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FIG. 1. Neutrino energy spectra for i/ M , i> M , and ^ e 
at the CCFR detector for the Fermilab wide-band neutrino 
beam. 

2. v e ,nN — > ^e iM X (f e ,^ neutral- current (NC) events). 

3. v e N — > eA (f e CC events). 

The majority (97.7%) of events observed in the de- 
tector are produced by muon neutrino interactions. The 
Vp CC events can be identified by the presence of a muon 
in the final state which penetrates beyond the end of 
the hadron shower, depositing energy characteristic of a 
minimum ionizing particle [p| in a large number of con- 
secutive scintillation counters. Conversely, the electron 
produced in a v e CC event deposits energy in a few coun- 
ters immediately downstream of the interaction vertex 
and is typically much shorter than the hadron shower. 
The separation of v e4l NC from the v e CC events is ac- 
complished by using the difference in energy deposition 
pattern within the shower region; the v e CC events have 
a larger fraction of their energy deposited near the shower 
vertex. 

In this analysis, the three most important experimental 
quantities calculated for each event are length, visible en- 
ergy, and shower energy deposition profile. Event length 
is determined to be the number of scintillation counters 
spanned from the event vertex to the last counter with 
greater than a minimum-ionizing pulse height. The visi- 
ble energy in the calorimeter, E vis , is obtained by sum- 
ming the energy deposited in scintillation counters from 
the interaction vertex to five counters beyond the end of 
the shower. The shower energy deposition profile is char- 
acterized by the ratio of the sum of the energy deposited 
in the first three scintillation counters to the total visible 
energy. Accordingly, we define 



i 



E\ + E2 + E3 



(2) 



where £7, is the energy deposited in the i th scintillation 
counter downstream of the interaction vertex. 

The event length is determined by the end of the 
hadron shower for NC and v e CC events but is de- 
termined by the muon track for most CC events. To 
isolate events without a muon track we parameterize the 
event length as a function of energy for which 99% of 
hadron showers are contained as 



4.0 + 3.81 xlog(£ ms ). 



(3) 



Events which deposit energy over an interval less than 
Ljvc counters are classified as "short", otherwise they 
are "long". The long event sample consists almost ex- 
clusively of class 1 events, while the short sample is a 
mixture of class 2, class 3, and class 1 events with a low 
energy muon. 

Events were selected with at least 30 GeV deposited 
in the target calorimeter to ensure complete trigger effi- 
ciency. Additionally, we require the event vertex to be at 
least five counters from the upstream end and more than 
£jvc + 5 counters from the downstream end of the target, 
and less than 127 cm from the detector center-line. The 
resulting data sample consists of 632338 long events and 
291354 short events. 

To directly compare the long and short events a muon 
track from the data was added to the short events to 
compensate for the absence of a muon in NC events. The 
fraction, /, of CC events with a low energy muon con- 
tained in the short sample (which will now contain two 
muon tracks) was determined from a Monte Carlo sim- 
ulation to be approximately 20%. A simulated sample 
of such events was obtained by choosing a long event 
with the appropriate energy distribution from the data 
and combining it with a second short muon track. The 
length of the short track and its angular distribution were 
obtained from a Monte Carlo sample of CC events. 

A sample of v e CC interactions with a muon track 
added were obtained by convolving an electromagnetic 
shower generated using GEANT with an event from 
the long data sample with the appropriate energy. This 
assumes — v e universality. The energy distribution of 
v^s and the fractional energy transfer y were obtained 
from Monte Carlo. Because the hadron showers in the 
long data sample already have a muon track, the v e CC 
sample can be compared directly with the short and long 
events. 

The long and short 773 distributions were further cor- 
rected by subtracting contamination due to cosmic ray 
events. The cosmic ray background was estimated from 
an event sample collected during a beam-off gate using an 
analysis procedure identical to the one used for the data 
gates. Additionally, the 773 distribution of short CC 
events, normalized to the predicted fraction /, was sub- 
tracted from the short event sample. The 7/3 distributions 
for short, long, and v e CC events for various energy bins 
are shown in Figure ||. 

We extract the number of v e CC events in each of 15 
E vis bins by fitting the corrected shape of the observed 773 
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FIG. 2. ?73 distributions for short (solid line), long (dashed 
line), and v £ CC (dotted line) events in four of the energy bins 
studied. The v £ CC and long distributions are normalized to 
the respective number of events predicted by the fit. 



distribution for the short sample to a linear combination 
of long and v e CC 773 distributions: 



r] 3 {short) = a rj 3 (long) + (3 r] 3 (u e CC) 



(4) 



The x 2 °f the fit in each of the 15 E V i S bins ranges from 
33 to 78 for 41 degrees of freedom (DoF) with a mean 
value of 48. 

To search for v e oscillations the measured absolute flux 
of VeS at the detector was compared to the flux pre- 
dicted by a detailed beamline simulation 14 1. Figure || 
shows the measured number of v e CC's for each energy 
bin compared with the predicted flux. The x 2 value with 
a no-oscillations assumption is 6.8/15 DoF. We interpret 
a deficit in the measured v e flux as v e — > v T (or v e — > v s ) 
oscillations. (y e — * oscillations are excluded above 
mixing angles of 2 x 10~ 3 at the 90% confidence level in 
this Am 2 range). 

If v e — > v T oscillation occurs, some fraction, f e , of v r 
charged-current interactions will be observed in our v e 
data sample. These v e CC-like events result from v T in- 
teractions in which a large fraction of energy deposited 
by the final state r is electromagnetic. To determine 
this fraction we simulated charged-current v r interac- 
tions in our detector using GEANT and a combination 
of LUND |ll| to generate charged-current neutrino inter- 
actions and TAUOLA JlJ] to simulate tau lepton decays. 
We fit the resulting v T charged-current Monte Carlo sam- 
ple to a linear combination of pure v e CC and v e ^ NC 
generated samples. The resulting v e CC-like fraction of 
v T CC events is 18% for our data sample. 
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FIG. 3. Number of electron neutrinos as a function of vis- 
ible energy. For electron neutrinos the visible energy is equal 
to the total neutrino energy. The filled band shows Monte 
Carlo prediction assuming no oscillations. The curves shown 
are the effect of v a — > v T oscillations for sin 2 2a = 1 and 
Am 2 = 150 eV 2 (dashed) and Am 2 = 10000 eV 2 (dotted). 



The effect of v e — > v T oscillations on the observed v & 
spectrum was determined in the following way: a beam- 
line simulation was used to tag the creation point of a 
v e along the decay pipe giving the survival probablity for 
each v e as (1 - P(v e -+v T )) from Eq. (§). The predicted 
v e flux was normalized to the observed charged-current 
muon neutrino flux at the detector which was simulated 
in the same beamline Monte Carlo. We also added in the 
number of v T charged-current interactions which would 
appear in the extracted v e sample at the detector. The 
probablity of observing a v T charged-current interaction 
in this data sample was calculated from the predicted 
normalized v e flux multiplied by the v T creation prob- 
ability, P{v e — > v T ), and the v e CC-like fraction, f e . 
We took into account effects of v T charged-current cross 
section suppression by including mass suppression terms 
p3| , kinematic suppression for massive particle produc- 
tion, and the altered visible energy spectrum (E vis ) for 
v T charged-current events which contains visible energy 
from the tau decay in determining the effect of v r ap- 
pearance. 



The effect of v„ 



oscillations on the v e spec- 



trum depends on the creation probability of v T from v^, 
P{ y ij, ~ * v t)- Tli e v T appearance effect is calculated by 
multiplying Piy^ — > v T ) by the v e CC-like fraction, f e 
and weighting by the v r CC cross section suppression 
factor. For completeness, we include a limit on — * v r 
from this data sample (see Figure ||). 

The major sources of uncertainties in the comparison 
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TABLE I. The result for sin 2 2a from the fit at each Am 2 
for v e — > v r oscillations. The 90% C.L. upper limit is equal 
to the best fit sin 2 2a + 1.28cr. 



m 2 (eV 2 ) 


Best fit 


a 


10.0 


-0.775 


2.046 


20.0 


-0.208 


0.553 


30.0 


-0.112 


0.269 


40.0 


-0.046 


0.180 


50.0 


-0.045 


0.134 


60.0 


-0.026 


0.109 


70.0 


-0.012 


0.096 


80.0 


-0.006 


0.086 


90.0 


0.003 


0.079 


100.0 


0.006 


0.078 


125.0 


0.005 


0.073 


150.0 


0.013 


0.072 


175.0 


0.029 


0.072 


200.0 


0.033 


0.071 


225.0 


0.047 


0.074 


250.0 


0.070 


0.080 



Am 2 (o\r 2 \ 
/ lev j 


R DC f flf 
J3CSL 11 L 




(7 


275.0 


0.083 





087 


300.0 


0.094 





089 


350.0 


0.065 





092 


400.0 


0.025 





098 


450.0 


-0.011 





104 


500.0 


-0.024 





111 


600.0 


0.045 





110 


700.0 


0.065 





112 


800.0 


0.045 





115 


1000.0 


0.040 





124 


1500.0 


0.045 





122 


2000.0 


0.062 





124 


5000.0 


0.053 





125 


10000.0 


0.060 





124 


20000.0 


0.048 





123 



of the v e flux extracted from the data to that predicted 
by the Monte Carlo are: the statistical error from the 
fit in v e flux extraction, error in shower shape model- 
ing (described below), uncertainty in the absolute en- 
ergy calibration of the detector (1%) which affects the 
relative neutrino flux extracted using a data sample with 
low hadron energy [fi5[ , and finally the uncertainty in the 
predicted flux of v e 's at the detector which is estimated 
to be 4.1% Q. This error is dominated by a 20% pro- 
duction uncertainty in the Kl content of the secondary 
beam which produces 16% of the v e flux. The majority 
of the v e flux comes from Kf^ decays, which are well- 
constrained by the observed spectrum from K^ 2 de- 
cays Other sources of systematic errors were also 
investigated and found to be small. 

The uncertainty in shower shape modeling is estimated 
by extracting the v e flux using two definitions of 77. Anal- 
ogous to the definition of 773 given in Eq. (^) , we define 1J4 
to be the ratio of the sum of the energy deposited outside 
the first four scintillation counters to the total visible en- 
ergy. If the modeling of the showers were correct, the 
difference in the number of electron neutrinos measured 
by the two methods should be small, any difference is 
used to estimate the systematic error. Since this error 
was shown not to be correlated among energy bins, we 
add it in quadrature to the statistical error from the fit. 

The data are fit by forming a x 2 which incorporates 
the Monte Carlo generated effect of oscillations and terms 
with coefficients accounting for systematic uncertainties. 
A best fit sin 2 2a is determined for each Am 2 by minimiz- 
ing the x 2 as a function of sin 2 2a and these systematic 
coefficients. At all Am 2 , the data are consistent with no 
observed oscillations. Table | shows the best fit value of 
sin 2 2a at each Am 2 for v e — > v T oscillations. The largest 




FIG. 4. Excluded region of sin 2a and Am for v £ — > v T 
and Ve — > u s oscillations from this analysis at 90% confidence 
(one-sided limit). 

statistical significance of a best-fit oscillation at any Am 2 
is la. 

The frequentist approach |l(| is used to set a 90% con- 
fidence upper limit for each Am 2 . The limit in sin 2 2a 



corresponds to a shift of 1.64 units in 
mum x 2 (at the best fit value in Table Jl 
fidence upper limit is plotted in Figure 



2 from the mini- 
1. The 90% con- 
fer v e — > v T and 



v e — > v s . The best limits of sin 2 2a are < 9.9 x 10 2 is at 
Am 2 = 125 eV 2 and < 8.3 x 10~ 2 is at Am 2 = 125 eV 2 
respectively. For sin 2 2a — 1, Am 2 > 20 eV 2 is excluded, 
and sin 2 2a > 0.21 for Am 2 > 1000 eV 2 for v e -> v r . 

As an alternative statistical treatment of this result we 
present 90% confidence limits based on the unified ap- 

17l] recently adopted by 



CCFR limits obtained 



proach of Feldman and Cousins 
the PDG U). Figure I shows al" 
using the longitudinal shower-shape method. Our pre- 
viously published limit [Q on — > v e used a one-sided 
confidence limit approach (as above). 

In conclusion, we have used a high-statistics sample 
of v e charged-current interactions in the CCFR coarse- 
grained calorimetric detector to search for v e — > v T , 



v s and v u 



oscillations. We see a result consis- 



tent with no neutrino oscillations and find 90% confidence 
level excluded regions in sin 2 2a — Am 2 phase space. This 



4 



result improves on existing limits for v e — > v T in the range 
50 eV 2 < Am 2 < 200 eV 2 . 
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sin 2a 



FIG. 5. Excluded region of sin 2 2a and Am 2 for (right 
to left) v e — * v T , — * v a (disappearance), i/ M — > v T , and 
— > v e (see Reference [7]) at 90% confidence using the Feld- 
man-Cousins approach. The first three limits are new. The 
v^ — > v £ limit differs from Reference [7] only in the construc- 
tion of the 90% confidence limit. 
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